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Abstract 

We consider dimer models on graphs which are bipartite, periodic and satisfy 
a geometric condition called isoradiality, defined in [18]. We show that the scaling 
limit of the height function of any such dimer model is 1/-ydr times a Gaussian free 
field. Triangular quadri-tilings were introduced in [6]; they are dimer models on 
a family of isoradial graphs arising form rhombus tilings. By means of two height 
functions, they can be interpreted as random interfaces in dimension 2 + 2. We 
show that the scaling limit of each of the two height functions is l/y^ times a 
Gaussian free field, and that the two Gaussian free fields are independent. 


1 Introduction 

1.1 Height fluctuations for isoradial dimer models 

1.1.1 Dimer models 

The setting for this paper is the dimer model. It is a statistical mechanics model repre¬ 
senting diatomic molecules adsorbed on the surface of a crystal. An interesting feature 
of the dimer model is that it is one of the very few statistical mechanics models where 
exact and explicit results can be obtained, see [14, 15] for an overview. Another very 
interesting aspect is the alleged conformal invariance of its scaling limit, which is al¬ 
ready proved in the domino and 60°-rhombus cases [16, 17, 19]. Theorem 1 of this 
paper shows this property for a wide class of dimer models containing the above two 
cases. 

In order to give some insight, let us precisely define the setting. The dimer model is in 
bijection with a mathematical model called the 2-tiling model representing discrete 
random interfaces. The system considered for a 2-tiling model is a planar graph G. 
Configurations of the system, or 2-tilings, are coverings of G with polygons consisting 
of pairs of edge-adjacent faces of G, also called 2-tiles, which leave no hole and don’t 
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overlap. The system of the corresponding dimer model is the dual graph G* of G. 
Configurations of the dimer model are perfect matchings of G*, that is set of edges 
covering every vertex exactly once. Perfect matchings of G* determine 2-tilings of G as 
explained by the following correspondence. Denote by /* the dual vertex of a face / of 
G , and consider an edge f*g* of G*. We say that the 2-tile of G made of the adjacent 
faces / and g is the 2-tile corresponding to the edge f*g*. Then 2-tiles corresponding 
to edges of a dimer configuration form a 2-tiling of G. Let us denote by A4(G*) the set 
of dimer configurations of G*. 


As for all statistical mechanics models, dimer configurations are chosen with respect to 
the Boltzmann measure defined as follows. Suppose that the graph G* is finite, and 
that a positive weight function n is assigned to edges of G* , then each dimer configura¬ 
tion M has an energy , £(M) = — ^ e eMl°§K e )- The probability of occurrence of the 
dimer configuration M chosen with respect to the Boltzmann measure /i 1 is: 


/(M) 


e -£(M) 

Z(G*,u) 


n c gM Ke) 

Z{G*,u) ’ 


where Z(G*,u) is the normalizing constant called the partition function. Using the 
bijection between dimer configurations and 2-tilings, v can be seen as weighting 2-tiles, 
and /r 1 as a measure on 2-tilings of G. When the graph G* is infinite, a Gibbs measure 
is defined to be a probability measure on A4 (G*) with the following property: if the 
matching in an annular region is fixed, then matchings inside and outside of the annulus 
are independent, moreover the probability of any interior matching M is proportional 
to n e eM u ( e )- F rom now on, let us assume that the graph G satisfies condition (*) 
below: 


(*) 


The graph G is infinite, planar, and simple ( G has no loops and no 
multiple edges); its vertices are of degree >3. G is simply connected , 
i.e. it is the one-skeleton of a simply connected union of faces; and it is 
made of finitely many different faces, up to isometry. 


1.1.2 Isoradial dimer models 

This paper actually proves conformal invariance of the scaling limit for a sub-family of 
all dimer models called isoradial dimer models, introduced by Kenyon in [18]. Much 
attention has lately been given to isoradial dimer models because of a surprising feature: 
many statistical mechanics quantities can be computed in terms of the local geometry 
of the graph. This fact was conjectured in [18], and proved in [7]. The motivation 
for their study is further enhanced by the fact that the yet classical domino and 60°- 
rhombus tiling models are examples of isoradial dimer models. Last but not least their 
understanding allows us to apprehend a random interface model in dimension 2 + 2 
called the triangular quadri-tiling model introduced in [6], see Section 1.2.1. 

Let us now define isoradial dimer models. Speaking in the terminology of 2-tilings, 
isoradial 2-tiling models are defined on graphs G satisfying a geometric condition 
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called isoradiality: all faces of an isoradial graph are inscribable in a circle, and 
all circumcircles have the same radius, moreover all circumcenters of the faces are 
contained in the closure of the faces. The energy of configurations is determined by a 
specific weight function called the critical weight function, see Section 2.1 for definition. 
Note that if G is an isoradial graph, an isoradial embedding of the dual graph G* is 
obtained by sending dual vertices to the center of the corresponding faces. Hence, the 
corresponding dimer model is called an isoradial dimer model. 


1.1.3 Height functions 

Let G be an isoradial graph, whose dual graph G* is bipartite. Then, by means of 
the height function, 2-tilings of G can be interpreted as random discrete 2-dimensional 
surfaces in a 3-dimensional space that are projected orthogonally to the plane. In 
physics terminology, one speaks of random interfaces in dimension 2 + 1. The height 
function, denoted by h, is an M-valued function on the vertices of every 2-tiling of G, 
and is defined in Section 3. 


1.1.4 Gaussian free field in the plane 

The Gaussian free field in the plane is defined in Section 4. It is a random distri¬ 
bution which assigns to functions £ C^M 2 ) (the set of compactly sup¬ 

ported smooth functions of M 2 , which have mean 0), a real Gaussian random vector 
(+Vi,..., Ftpk) whose covariance function is given by 


E(FipiFipj) = 


g(x,y)g>i(x)tpj(y)dx dy, 


where g(x,y) = —^log|x — y\ is the Green function of the plane (defined up to an 
additive constant). The Gaussian free field is conformally invariant [17]. 


1.1.5 Statement of result 

Let G be an isoradial graph, whose dual graph G* is bipartite. Suppose moreover that 
G* is doubly periodic, i.e. that the graph G* and its vertex-coloring are periodic. Then 
by Sheffield’s theorem [26], there exists a two-parameter family of translation invariant, 
ergodic Gibbs measures; let us denote by y the unique measure which has minimal free 
energy per fundamental domain. From now on, we assume that dimer configurations 
of G* are chosen with respect to the measure y. 

Let us multiply the edge-lengths of the graph G by e > 0, this yields a new graph G s . 
Let h £ be the unnormalized height function on 2-tilings of G £ . An important issue in 
the study of the dimer model is the understanding of the fluctuations of h £ , as the mesh 
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£ tends to 0. This question is answered by Theorem 1 below. Define: 

H £ : C^q(M 2 ) -»• K 

+ i—* H £ ip = e 2 a(v*)<p(v)h £ (v), 

veV{G q 

where V(G £ ) denotes the set of vertices of the graph G £ , and a(w*) is the area in G* of 
the dual face v* of a vertex v. 

Theorem 1 Consider a graph G satisfying the above assumptions, then H £ converges 
weakly in distribution to -A times a Gaussian free field, that is for every <pi ,..., £ 

CT^(M 2 ), (H £ p i,... ,H £ ipk) converges in law (as e —> 0 ) to -^(Ftpi,... ,Fp>h), where 
F is a Gaussian free field. 

• As a direct consequence of Theorem 1, we obtain convergence of the height function 
of domino and 60°-rhombus tilings chosen with respect to the uniform measure to a 
Gaussian free field. Note that this result is slightly different than those of [16, 17, 19] 
since we work on the whole plane, and not on simply connected regions. 

• The method for proving Theorem 1 is essentially that of [16], except Lemma 20 which 
is new. Nevertheless, since we work with a general isoradial graph (and not the square 
lattice), each step is adapted in a non trivial way. 

1.2 The case of triangular quadri-tilings 
1.2.1 Triangular quadri-tiling model 

An exciting consequence of Theorem 1 is that it allows us to understand height fluctu¬ 
ations in the case of a random interface model in dimension 2 + 2, called the triangular 
quadri-tiling model. It is the first time this type of result can be obtained on such a 
model. 

Let us start by defining triangular quadri-tilings. Consider the set of right triangles 
whose hypotenuses have length 1, and whose interior angle is n/3. Color the vertex at 
the right angle black, and the other two vertices white. A quadri-tile is a quadrilateral 
obtained from two such triangles in two different ways: either glue them along a leg 
of the same length matching the black (white) vertex to the black (white) one, or glue 
them along the hypotenuse. There are four types of quadri-tiles classified as I, II, III, 
IV, each of which has four vertices, see Figure 1 (left). A triangular quadri-tiling of 
the plane is an edge-to-edge tiling of the plane by quadri-tiles that respects the coloring 
of the vertices, see T of Figure 1 for an example. Let Q denote the set of all triangular 
quadri-tilings of the plane, up to isometry. 

In [6], triangular quadri-tilings of Q are shown to correspond to two superposed dimer 
models in the following way, see also Figure 1. Define a lozenge to be a 60°-rhombus. 
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Then triangular quadri-tilings are 2-tilings of a family of graphs £ which are lozenge- 
with-diagonals tilings of the plane, up to isometry. Indeed let T 6 Q be a triangular 
quadri-tiling, then on every quadri-tile of T draw the edge separating the two right tri¬ 
angles, this yields a lozenge-with-diagonals tiling L(T) called the underlying tiling. 
Moreover, the lozenge tiling L(T) obtained from L[T ) by removing the diagonals, is a 
2-tiling of the equilateral triangular lattice T. 


T is a 2-tiling of L (T) 



L(T) is a 2-tiling of the triangular lattice 


Figure 1: Four type of quadri-tiles (left). Triangular quadri-tilings correspond to two 
superposed dimer models (right). 

Note that lozenge-with-diagonals tilings and the equilateral triangular lattice are iso- 
radial graphs. Assigning the critical weight function to edges of their dual graphs, we 
deduce that triangular quadri-tilings of Q correspond to two superposed isoradial dimer 
models. 

1.2.2 Height functions for triangular quadri-tilings 

Triangular quadri-tilings are characterized by two height functions in the following way. 
Let T £ Q be a triangular quadri-tiling, then the first height function, denoted by hi, 
assigns to vertices of T the “height” of T interpreted as a 2-tiling of its underlying 
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lozenge-with-diagonals tiling L(T). The second height function, denoted by / 12 , assigns 
to vertices of T the height of L(T) interpreted as a 2-tiling of T. An example of 
computation is given in Section 3.3. By means of h\ and h 2 , triangular quadri-tilings 
are interpreted in [6] as discrete random 2-dimensional surfaces in a 4-dimensional space 
that are projected orthogonally to the plane, i.e. in physics terminology, as random 
interfaces in dimension 2 + 2. 

1.2.3 Statement of result 

The notion of Gibbs measure can be extended naturally to the set Q of all triangular 
quadri-tilings, see Section 2.4. In [7], we give an explicit expression for such a Gibbs 
measure P, and conjecture it to be of minimal free energy per fundamental domain 
among a four -parameter family of translation invariant, ergodic Gibbs measures. Let 
us assume that triangular quadri-tilings of Q are chosen with respect to the measure 
P. 

Corollary 2 below describes the fluctuations of the unnormalized height functions h\ 
and /if. Suppose that the equilateral triangular lattice has edge-lengths 1, and let T £ 
be the lattice T whose edge-lengths have been multiplied by e. Observe that vertices 
of T are vertices of L, for every lozenge-with-diagonals tiling L E C. For i = 1,2, and 
for ip € (M 2 ) define: 


Hfp = e 2 i( v )> 

v£V(T s ) 

Corollary 2 For i = 1, 2, and every p\, . .. ,pk € C'^(M 2 ), (Hfcpi ,..., Hfp &) con¬ 
verges in law (as e —► 0) to -^(Fnpi,... ,Fi(pk), where F\ is a Gaussian free field. 
Moreover, F\ and Ff are independent. 

1.3 Outline of the paper 

• Section 2: statement of the explicit expressions of [6] for the Gibbs measures p 
and P, that are used in the proof of Theorem 1 and Corollary 2. 

• Section 3: definition of the height function on vertices of 2-tilings of isoradial 
graphs. 

• Section 4: definition of the Gaussian free field of the plane. 

• Section 5 and Section 6: proof of Theorem 1 and Corollary 2. 

Acknowledgments: We would like to thank Richard Kenyon for proposing the questions 
solved in this paper, and for the many enlightening discussions. We are grateful to 
Erwin Bolthausen, Cedric Boutillier and Wendelin Werner for their advice and sugges¬ 
tions. 
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2 Minimal free energy Gibbs measure for isoradial dimer 
models 

In the whole of this section, we let G be an isoradial graph, whose dual graph G* is 
bipartite; B denotes the set of black vertices, and W the set of white ones. In the proof 
of Theorem 1, we use the explicit expression of [6] for the minimal free energy per 
fundamental domain Gibbs measure n on 2-tilings of G, and in the proof of Corollary 
2, we use the explicit expression of [6] for the Gibbs measure P on triangular quadri- 
tilings. The goal of this section is to state the expressions for /j and P. In order to 
do so, we first define the critical weight function and the Dirac operator, introduced in 
[18]- 

2.1 Critical weight function 

2.1.1 Definition 

The following definition is taken from [18]. To each edge e of G*, we associate a unit 
side-length rhombus R(e) whose vertices are the vertices of e and of its dual edge e* 
(R(e) may be degenerate). Let R = U ee G*f?(e). Tbe critical weight function v at 
the edge e is defined by v(e) = 2 sin 8, where 28 is the angle of the rhombus R{e) at the 
vertex it has in common with e; 8 is called the rhombus angle of the edge e. Note 
that v(e) is the length of e*. 


2.1.2 Example: critical weights for triangular quadri-tilings 

Recall that triangular quadri-tilings of Q correspond to two superposed isoradial dimer 
models, the first on lozenge-with-diagonals tilings and the second on the equilateral 
triangular lattice T. Let us note that the dual graphs of lozenge-with-diagonals tilings 
and of T are bipartite. We now compute the critical weights in the above two cases. 

Consider the equilateral triangular lattice T, then edges of its dual graph T*. known as 
the honeycomb lattice, all have the same rhombus angle, equal to 7r/3, and the same 
critical weight, equal to \/3. 

Consider a lozenge-with-diagonals tiling L 6 C. Observe that the circumcenters of the 
faces of L are on the boundary of the faces, so that in the isoradial embedding of the 
dual graph L* some edges have length 0, and the rhombi associated to these edges are 
degenerate, see Figure 2. Since edges of L* correspond to quadri-tiles of L , we classify 
them as being of type I, II, III, IV. Figure 2 below gives the rhombus angles and the 
critical weights associated to edges of type I, II, III and IV, denoted by ei, en, em, eiv 
respectively. 
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Figure 2: Dual graph of a lozenge-with-diagonals tiling (left). Critical weights for 
quadri-tiles (right). 

2.2 Dirac and inverse Dirac operator 

Results in this section are due to Kenyon [18], see also Mercat [25]. Define the Hermitian 
matrix K indexed by vertices of G* as follows. If v\ and V 2 are not adjacent, K(y 1 ,^ 2 ) = 
0. If w £ IV and b G B are adjacent vertices, then K(w,b) = K{b,w) is the complex 
number of modulus v(wb) and direction pointing from w to b. Another useful way to 
say this is as follows. Let R(wb) be the rhombus associated to the edge wb, and denote 
by w,x,b,y its vertices in cclw (counterclockwise) order, then K(w,b ) is i times the 
complex vector x — y. If w and b have the same image in the plane, then | I\(w, 6)| = 2, 
and the direction of K(w,b ) is that which is perpendicular to the corresponding dual 
edge, and has sign determined by the local orientation. The infinite matrix K defines 
the Dirac operator K: by 

( Kf)(v)= Y K i v , u )f( u ), 

u£G* 

where V(G*) denotes the set of vertices of the graph G*. 

The inverse Dirac operator K~ 1 is defined to be the operator which satisfies 

1. KK~ l = Id, 

2. JL _1 (6, w) —> 0, when \b — w\ —> 00 . 

In [18], Kenyon proves uniqueness of and existence by giving an explicit expression 
for K~ l {b,w) as a function of the local geometry of the graph. 

2.3 Minimal free energy Gibbs measure for isoradial graphs 

If e\ = w\b\,... , e*, = Wkbk is a subset of edges of G*, define the cylinder set 
{ei,... ,ek} of G* to be the set of dimer configurations of G* which contain the edges 
ei,..., e#;. Let A be the field consisting of the empty set and of the finite disjoint unions 
of cylinders. Denote by cr(A) the cr-field generated by A. 
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Theorem 3 [6] Assume G* is doubly periodic. Then, there is a probability measure p 
on (. M.{G*),a(A )) such that for every cylinder {ei,..., e^} of G*, 


h( e li ■■■■> e k) 



(1) 


Moreover p is a Gibbs measure on M(G*), and it is the unique Gibbs measure which 
has minimal free energy per fundamental domain among the two-parameter family of 
translation invariant, ergodic Gibbs measures of [26]. 


Remark 4 

• Refer to [21] for the definition of the free energy per fundamental domain. 

• In [6], we prove that the periodicity assumption can be released in the case of lozenge- 
with-diagonals tilings. That is, given any lozenge-with-diagonals tiling of £, equation 
(1) defines a Gibbs measure on dimer configurations of its dual graph. Although funda¬ 
mental domains make no sense in case of non-periodic graphs, the minimal free energy 
property can still be interpreted in some wider sense. 


2.4 Gibbs measure on triangular quadri-tilings 

The construction of this section is taken from [6]. Consider the set Q of all triangular 
quadri-tilings of the plane up to isometry, and assume that quadri-tiles are assigned 
a positive weight function. Then the notion of Gibbs measure on Q is a natural 
extension of the one used in the case of dimer configurations of fixed graphs. It is a 
probability measure that satisfies the following: if a triangular quadri-tiling is fixed 
in an annular region, then triangular quadri-tilings inside and outside of the annulus 
are independent; moreover, the probability of any interior triangular quadri-tiling is 
proportional to the product of the weights of the quadri-tiles. Denoting by JA the set 
of dimer configurations corresponding to triangular quadri-tilings of Q, and using the 
bijection between Q and A4, we obtain the definition of a Gibbs measure on A4. 

Define £* to be set of dual graphs L* of lozenge-with-diagonals tilings L G £. Although 
some edges of £* have length 0, we think of them as edges of the one skeleton of the 
graphs, so that to every edge of £*, there corresponds a unique quadri-tile. Let e be an 
edge of £*, and let q e be the corresponding quadri-tile, then q e is made of two adjacent 
right triangles. If the two triangles share the hypotenuse edge, they belong to two 
adjacent lozenges; else if they share a leg, they belong to the same lozenge. Let us 
call these lozenge(s) the lozenge(s) associated to the edge e, and denote it/them 
by l e (that is l e consists of either one or two lozenges). Let k e be the edge(s) of T* 
corresponding to the lozenge(s) l e . Let us introduce one more definition, if {e\, ..., e*,} 
is a subset of edges of £*, then the cylinder set {ei,...,efc} is the set of dimer 
configurations of A4 which contain these edges. Denote by C the field consisting of 
the empty set and of the finite disjoint unions of cylinders. Denote by <r(C) the cr-field 
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generated by C. 

Consider a lozenge-with-diagonals tiling L £ C, and denote by /u L the minimal free 
energy per fundamental domain Gibbs measure on (J\4(L*),a(A)) given by Theorem 
3, where a (A) is the u-field of cylinders of Similarly, denote by /r T the minimal 

free energy per fundamental domain Gibbs measure on (At(T*), ct(£>)), where a(B) is 
the a-field of cylinders of At(T*). 

Let us define fi L on (A i,a(C)) by: 

[ H L {e i, ...,e k ) if the lozenges l ei ,..., l efc belong to L, 

A (ei, • • • ,e k ) = < 

[ 0 else, 

where we recall that L is the lozenge tiling obtained from the lozenge-with-diagonals 
tiling L by removing the diagonals. Then, it is easy to check that fi L is a probability 
measure on (Af, o(C)). In order to simplify notations, we write /u L for fi L whenever no 
confusion occurs. 

Now, on (Af x Af(T*),£> x C), define: 

P((ei,... ,e k ) x (ki, ..., k m )) = ^ ^ L (ei,..., e k )dn T (l*). 

{L*eX(T*):fci,...,fc m eL*} 

Using Kolomogorov’s extension theorem, P extends to a probability measure on (Af x 
T *,a{B x C)). Let us also denote by P the marginal of P on Af, then in [6], P is 
shown to be a Gibbs measure on Af, and conjectured to be of minimal free energy per 
fundamental domain among a four- parameter family of translation invariant, ergodic 
Gibbs measures. 

3 Height functions 

In the whole of this section, we let G be an isoradial graph whose dual graph G* 
is bipartite; as before, B denotes the set of black vertices, W the set of white ones. 
We define the height function h on vertices of 2-tilings of G , whose fluctuations are 
described in Theorem 1. As in [21], see also [4], h is defined using flows. 

The bipartite coloring of the vertices of G* induces an orientation of the faces of G: 
color the dual faces of the black (white) vertices black (white); orient the boundary 
edges of the black faces cclw, the boundary edges of the white faces are then oriented 
cw. 


3.1 Definition 

Let us first define a flow ojq on the edges of G*. Consider an edge wb of G*, then R(wb) 
is the rhombus associated to wb, and 9 w b is the corresponding rhombus angle. Define 
loq to be the white-to-black flow, which flows by d^/ir along every edge wb of G*. 
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Lemma 5 The flow ujq has divergence 1 at every white vertex, and — 1 at every black 
vertex of G*. 

Proof: 

By definition of the rhombus angle, we have 

Vw G W, ^ 29 w b = 27r; V 6 G B, ^ 29 w b = 27r. 

b:b~w wr.w^b 

□ 

Now, consider a 2-tiling T of G, and let M be the corresponding perfect matching of 
G*. Then M defines a white-to-black unit flow ui on the edges of G*: flow by 1 along 
every edge of M , from the white vertex to the black one. The difference uiq — to is a 
divergence free flow, which means that the quantity of flow that enters any vertex of 
G* equals the quantity of flow which exists that same vertex. 

We are ready for the definition of the height function h. Choose a vertex vq of G, 
and fix h(v o) = 0. For every other vertex v of T, take an edge-path 7 of G from vo to 
v. If an edge uv of 7 is oriented in the direction of the path, and if we denote by e its 
dual edge, then h increases by LOo(e) — ui{e) along uw, if an edge uv is oriented in the 
opposite direction, then h decreases by the same quantity along uv. As a consequence 
of the fact that wo - w is a divergence free flow, the height function h is well defined. 

The following lemma gives a correspondence between height functions defined on ver¬ 
tices of G , and 2-tilings of G. 

Lemma 6 Let h be an M-valued function on vertices of G satisfying 

• h(v 0 ) = 0 , 

• h(y) — h(u ) = iOo(e) or u;o(e) — 1 for any edge uv oriented from u to v, where e 
denotes the dual edge ofuv. 

Then, there is a bijection between functions h satisfying these two conditions, and 2- 
tilings of G. 

Proof: 

The idea of the proof closely follows [8]. Let T be a 2-tiling of G , M be the corresponding 
matching, and w be the unit white-to-black flow defined by M. Then, the height 
function h satisfies the conditions of the lemma: consider an edge uv of G oriented 
from u to v and denote by e its dual edge, then h{v) — h(u ) = wo(e) — w(e), and by 
definition cj(e) = 0 or 1 . 

Conversely, consider an M-valued function h as in the lemma. Let us construct a 2 -tiling 
T whose height function is h. Consider a black face F of G, and let e \,..., e m be the 
dual edges of its boundary edges. Then X^i w o(ei) = 1, so that there is exactly one 
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boundary edge uv along which h(v) — h(u) is ujq (e*) — 1 (where e* is the dual edge of uv). 
To the face F. we associate the 2-tile of G which is crossed by the edge uv. Repeating 
this procedure for all black faces of G, we obtain T. □ 

3.2 Interpretation 

The discrete interface interpretation of 2-tilings was first given by Thurston in the case 
of lozenges [30]. Following him, a 2-tiling of G can be seen as a discrete 2-dimensional 
surface S in a 3-dimensional space that has been projected orthogonally to the plane; 
the “height” of S is given by the function h. Stated in physics terminology, 2-tilings of 
G are random interfaces in dimension 2 + 1. 

3.3 Example: the two height functions of triangular quadri-tilings 

Consider a triangular quadri-tiling T £ Q. Then, recall that h\ assigns to vertices of 
T the “height” of T interpreted as a 2-tiling of its underlying lozenge-with-diagonals 
tiling L(T), and h -2 assigns to vertices of T the “height” of L(T) interpreted as a 2-tiling 
of T. Let us now explicitly compute h\ and h- 2 - 

The definition of the flow ujq(L(T)) on edges of L(T)* uses the rhombus angles of the 
edges of L(T)*. These have been computed in Section 2.1.2 and are equal to f, f, f 
for edges of type I, II, III and IV respectively. Hence if uv is the boundary edge of 
a quadri-tile of T, oriented from u to v, then the height change of h\ along uv is 
jjr, | depending on whether the dual edge e of the edge uv is of type I, II, III or IV 
respectively, see Figure 3. 

In a similar way, the flow woCT*) on edges of T* flows by | along every edge, so that 
if u'v' be a boundary edge of a lozenge of L(T) oriented from v! to v', then the height 
change of /12 along u'v' is equal to |. Thinking of a lozenge as the projection to the 
plane of the face of a cube [30], there is a natural way to assign a value for h ,2 at the 
vertex at the crossing of the diagonals of the lozenges, see Figure 3. 

By means of hi and h?, a triangular quadri-tiling T of Q is interpreted in [6] as a 
2-dimensional discrete surface S\ in a 4-dimensional space that has been projected 
orthogonally to the plane. >Sj can also be projected to ^Z 3 (Z 3 is the space Z 3 where 
the cubes are drawn with diagonals on their faces), and one obtains a surface S 2 ■ When 
projected to the plane S 2 is the underlying lozenge-with-diagonals tiling L{T). This 
can be restated by saying that triangular quadri-tilings of Q are discrete interfaces in 
dimension 2 + 2. 


12 





T 


L (T) 


0 

0 



T 


Figure 3: From left to right: triangular quadri-tiling T and first height function h\, 
underlying lozenge tiling L(T) and second height function h 2 , T with height functions 
hi (above) and I 12 (below). 


4 Gaussian free field of the plane 

Theorem 1 and Corollary 2 prove convergence of the height function h to a Gaussian 
free field. The goal of this section is to define the Gaussian free field of the plane. We 
refer to [9, 27] for other ways of the defining the Gaussian free field. 

4.1 The Green function of the plane, and Dirichlet energy 

The Green function of the plane, denoted by g, is the kernel of the Laplace equation 
in the plane, it satisfies A X g(x,y) = 6 x (y), where S x is the Dirac distribution at x. Up 
to an additive constant, g is given by 

9 (x,y) = ~^\og\x - y\. 

Define the following bilinear form 

G : C~ 0 (M 2 ) x C“ 0 (M 2 ) — ♦ R 

(¥> 1 ,^ 2 ) 1 —> G(ipi,ip 2 )= / g(x,y)(pi(x)ip 2 (y)dxdy. 

J R 2 J K 2 

G((p, (p) is called the Dirichlet energy of (p. Let us consider the topology induced by 
the L°° norm on (7^(M 2 ). 

Lemma 7 G is a continuous, positive definite, bilinear form. 

Proof: 

G is continuous 
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This is a consequence of the fact that for every pi,p 2 £ C)?g(M 2 ), the function g(x,y ) 
is integrable. 

G is positive definite 

For % = 1,2, denote by Kj = supp(^), and let ffix) = f R2 g{x,y)pfiy)dy. Let us prove 
that 

G(pi,p 2 ) = [ V/i(x).V/ 2 (x)dx. ( 2 ) 

J K 2 

For any R> 0, Green’s formula implies 

f Vffix).Vf 2 (x)dx = - [ Affix) f 2 (x)dx + f f 2 (x)V ffix).n(x)ds. (3) 

JB(0,R) Jb(0,R ) JS(0,R) 

Assume R is large enough so that K\ , K 2 C B(0,R). The first term of the right hand 
side of (3) satisfies 


/ Affix) ffix)dx = - A* ( / g(x,y)(pi(y)dy)[ g(x,y)p 2 (y)dy) dx, 

Jb(o,r) Jb(o,r) \Jk i / \Jk 2 / 


/ / g(x,y)pfix)p 2 {y)dydx = G{p l ,p 2 ). 


k i Jk 2 


In order to evaluate the second term of the right hand side of (3), let us compute 

V/i(x).n(x) = — / ¥ J i(y)Vlog|x-y|.n(x)dy, 

2tt Ja'j 

= ^ , 2 dy, 

2tt J Kl \x - y |- 

=- f pfiy)(-, — — in ~ i— dy - [ p\(y )-.— rdy, 

2ir J K / ny ’ \ \x-y\ 2 |x| J 2tt J K J ny> \x\ y ’ 

=- f oJi (y) ( t— - - | 0 — i— dy (since p\ is a mean 0 function). 

2 nJ Kl \\x-y I 2 \x\J 

Vx G S(0,R), V y G A'i, we have - |yy = O (^), hence | V/i(x).n(x)| = O (^j); 

Vx G S(0,R), we also have |/ 2 (x)| = OfiogR), thus the second term of the right hand 
side of (3) is O ■ Taking the limit as R —> oo in (3), we obtain (2). 


Let us assume G(p i, pfij = 0. By equality (2) this is equivalent to f R2 |V/i(x)| 2 cix = 0, 
hence Vfi = (0,0). Since pfix) = Affix) = div(V/i(x)), we deduce (pi = 0. □ 


4.2 Random distributions 

The following definitions are taken from [10]. A random function F associates to 
every function p G C^q(M 2 ) a real random variable Fp. For p\ ,..,, p^ G C^(M 2 ), 
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we suppose that the joint probabilities a n < Ftp n < b n , 1 < n < k are given, and 
we ask that they satisfy the compatibility relation. A random function F is linear if 

Vw,^eQ(K 2 ), 

F(anpi + (3tp 2 ) = aFip i + /3Fip 2 . 

It is continuous if convergence of the functions (p nj to ipj (1 < j < k ) implies 

lirn (Fip ni ,...,Fip n ) = (Fip 1 ,...,F(p k ), 

n—>oo 

that is, if P(x) (resp. P n (x)) is the probability measure corresponding to the random 
variable (Ftpi,... ,Ftp k ) (resp. (Fip ni ,... ,Ftp nk )), then for any bounded continuous 
function / 

lirn / f(xi,...,x k )dP n (x)= / f(xi,...,x k )d,P(x). 

n^ocj J 

A random distribution F is a random function which is linear and continuous. It is 
said to be Gaussian if for every linearly independent functions ip \,..., <p k € 
the random vector ( Fipi ,..., F<p k ) is Gaussian. 

Two random distributions F and G are said to be independent if for any functions 
<pi,...,ip k 6 C'^(M 2 ), the random vectors (Fipi ,..., Fip k ) and (Gipi ,..., Gtp k ) are 
independent. 


4.3 Gaussian free field of the plane 

Theorem 8 [3] If G : C(?jj(M 2 ) x C^M 2 ) —> M is a bilinear, continuous, positive 
definite form, then there exists a Gaussian random distribution F, whose covariance 
function is given by 

E(F(piFip 2 ) = G(tpi,ip 2 ). 


Using Lemma 7, and Theorem 8 , we define a Gaussian free field of the plane to be 
a Gaussian random distribution whose covariance function is 


E(F(piFtp 2 ) 


1 

2tt 


log \x - y\<pi(x)ip 2 {y) dx dy. 


5 Proof of Theorem 1 

We place ourselves in the context of Theorem 1: G is an isoradial graph whose dual 
graph G* is doubly periodic and bipartite; edges of G* are assigned the critical weight 
function, and G £ is the graph G whose edge-lengths have been multiplied by e. Recall 
the following notations: H e (p = e 2 Yl v eV(G s ) a ( v *) l P( v )h £ (v), /j is the minimal free en¬ 
ergy per fundamental domain Gibbs measure for the dimer model on G* of Section 2.3, 
and F is a Gaussian free held of the plane. 

Since the random vector (Fipi,... ,Fip k ) is Gaussian, to prove convergence of (H £ ip\, 
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... ,H £ (fk ) to (F<pi ,..., Fipk), it suffices to prove convergence of the moments of 
(H £ ipi,... ,H £ tpk) to those of (Fipi ,..., F<pk)\ that is we need to show that for ev¬ 
ery /c-tuple of positive integers (mi,..., m*,), we have 

limE[(fTVi) roi • • • (H £ <p k ) mk ] = E[(F^) mi ... (Fip k ) mk ). (4) 

£—>0 

In Section 5.1, we prove two properties of the height function h, and in Section 5.2, we 
give the asymptotic formula of [18] for the inverse Dirac operator K~ l . Using these 
results in Section 5.3, we prove a formula for the limit (as e —0) of the k th moment 
of h. This allows us to show convergence of E [{H £ (p) k \ to E[(T 1 <^)) fc ] in Section 5.4. One 
then obtains equation (4) by choosing ip to be a suitable linear combination of the <pi s. 

As before, B denotes the set of black vertices of G*, and W the set of white ones. 
Moreover, we suppose that faces of G have the orientation induced by the bipartite 
coloring of the vertices of G*. 

5.1 Properties of the height function 

Let u, v be two vertices of G , and let 7 be an edge-path of G from u to v. First, 
consider edges of 7 which are oriented in the direction of the path, that is edges which 
have a black face of G on the left, and denote by f\,, f n their dual edges. Hence 
an edge fj consists of a black vertex on the left of 7 , and of a white one on the right. 
Similarly, consider edges of 7 which are oriented in the opposite direction, and denote 
by ei,... ,e m their dual edges, hence an edge e* consists of a white vertex on the left 
of 7 , and of a black one on the right. Let I e be the indicator function of M(G*): 
I e (M) = 1, if the edge e belongs to the dimer configuration M of G*, and 0 else. 

Lemma 9 

m n 

h{v) - h{u) = ^(le, - Me?)) + + M/?))- 

3 = 1 4=1 


Proof: 

Let ej be the dual edge of an edge UjVj of 7 oriented from Vj to Uj. Denote by 9j the 
rhombus angle of the edge ej, then by Lemma 6 , 


h(yj ) — h(uj) 





7T 


if the edge ej belongs to the dimer configuration of G*, 
else. 


Hence h(yj) — h(uj ) = (— + l)I ej — ^-(1 — I ej ) = l ej — %■ Moreover, a direct 

computation using formula ( 1 ) yields Me.?) = 9j/ir, so 


h(vj) - h(uj) = I ej - n(ej). 
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Similarly, when fj is the dual edge of an edge u'jv'j of 7 oriented from u' 3 to u', we 
obtain h(vj ) — h(uj) = —If j + / J,(fj ), and we conclude 

ran m n 

h{v)-h{u) =^2h(v j )-h(u j )+J2 h {v , j )-h{u' j ) = +/r(/ i )). 

i=i i=i i=i i=i 

□ 


Lemma 10 

IE n[h{v) - h(u)] = 0. 

Proof: 

m n 

By Lemma 9 we have, E ^[h{v)-h{u)\ = ^E M [I ej .-/r(e :) -)]+^E jU [-I= 0. □ 

3 =1 i=i 


5.2 Asymptotics of the inverse Dirac operator A' 1 

In order to state the asymptotic formula of [18] for the inverse of the Dirac operator K 
indexed by vertices of G* , we let w £ W be a white vertex of G*, v any other vertex of 
G* , and define the rational function f wv (z) of [18]. Recall that R is the set of rhombi 
associated to edges of G* , and consider w = vo, v\ , V 2 , ■ ■ ■ , = v an edge-path of R 

from w to v. Each edge VjVj+i has exactly one vertex of G* (the other is a vertex of 
G). Direct the edge away from this vertex if it is white, and towards this vertex if it 
is black. Let e lolj be the corresponding vector in R (which may point contrary to the 
direction of the path). Then, f wv is defined inductively along the path, starting from 

fwU 0) = 1 . 


If the edge leads away from a white vertex, or towards a black vertex, then 

„ ( \ — f wv i (z) 

Jwvj+i [Z ) ^ g ioij ’ 

else, if it leads towards a white vertex, or away from a black vertex, then 

fwv j+ 1 ( 0 ) = fw Vj (z).(z - e iaj ). 

The function f wv (z ) is well defined (i.e. independent of the edge-path of R from w to 
v). Then, Kenyon gives the following asymptotics for the inverse Dirac operator K~ l \ 


Theorem 11 [18] Asymptotically, as \b — w\ —> 00 , 


K~\b,w) = 


1 + fwb( 0 ) 


1 


2tt V b — w b — w 


) +0 ( 


\b — u)| 
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5.3 Moment formula 

Let u\,..., Uk, v \,..., Vk be distinct points of M 2 , and let 71 ,, 7 ^ be pairwise disjoint 
paths such that 7 j runs from Uj to Vj. Let u £ j,Vj be vertices of G £ lying within O(e) of 
Uj and Vj respectively. In order to simplify notations, we write h for the unnormalized 
height function h £ on 2-tilings of G £ . Then, we have 

Proposition 12 For every k £ N, k > 2 


lim ¥.[(h(vl)-h{u\))... (h(v%) - h{u%))\ 
£ —>0 


where z® = Zi and z\ = Z{. 


(- 0 * 

( 27 r) fc 


E (- 1 )" 5 

5=0,1 


( 

[ . 

. / det 

( 5 ) dzf ... dz% 

J'Yl 

k i,je[i,k\ 

V Z i - Z j J 

V 


/ 


(5) 


Proof: 

Steps of the proof follow [16], but since we work in a much more general setting, they 
are adapted in a non-trivial way. 

Let 7f, • • • ,7^ be pairwise disjoint paths of G £ , such that 7 ? runs from u £ to v £ - and 

approximates 7 j within O(e). For every j, denote by fj s the dual edge of the s th edge 
of the path 7 ?, which is oriented in the direction of the path: fj s consists of a black 
vertex on the left of 7 1, and of a white one on the right. Denote by ejt the dual edge 
of the t th edge of the path 7 J, which is oriented in the opposite direction: e 3 t consists 
of a black vertex on the right of 7 1, and of a white one on the left. Using Lemma 9, we 
obtain 

E[(h(uf) - h{u\))... {h(v%) - h(u £ k ))] = 


= E 


E( Ie i*i -Meiti)) - -U/isi)) 


e kt k )) V{fks k )) 


— y . ^[Ie ltl M( e ltl )] • • • pkfct fc v( e kt k )] • ■ • + ( 1) y ] E[Iy lsi Ll(flsi )] . . . V-(fks 

s lvi s fe 


E E (-l) 5 l+ - + 5 ''E[Ie 4i -n{e tl )]...\ 




( 6 ) 


<5i,-..,<5fcG{0,l} t <5i t S k 


where t- j = * 3 J 


t-j if 5j = 0 


Sj if Sj = 1 


Uj 


J'V 


_ Ji _ J e jtj if ^ = 0 


fjsj if Sj — 1 


For the time being, let us drop the second subscript. Write ej = Wjbj and fj = w'jb'j. 
Moreover, let us introduce the notation Wj j , where Wj J = Wj if Sj = 0, and vj/ 1 = w' rj if 
Sj = 1, similarly we introduce the notation l/- J . Hence we can write a generic term of 


( 6 ) as 


( _ 1 ) *+...-w. E[(Im , lt? _ ... ( 1 ^ - ,*££))). 
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Lemma 13 [13] 


( _ 1) 5 1 +...+4 E[(I ^ ] 

b h ~v{w S ib S i)) 

■ ■ ■ (Vkft'k “ V( W k kb l k ))\ = 

k k 



0 

K~\b[\w^) 

K-\b{\w S *) 

= (—1 Y 1+ - +&k a E 

K~\b 5 2 \w[') 

0 

K- l {bt!,w S k “) 


K-\b 5 k \wl 1 ) 

k- 1 

0 

k 



(7) 

where a E = ]^[ K {w - 3 , b - ), and K 

is the Dirac operator indexed by vertices of G*. 


3 =1 


A typical term in the expansion of (7) is 

(~l) Sl+ - +Sk aE sgn aK~ l {b\ l ,w 5 J^)... K~ l (b { k , w S J$j ), ( 8 ) 

where a £ Sk, and Sk is the set of permutations of k elements, with no fixed points. 
To simplify notations, let us assume a is a fc-cycle, hence ( 8 ) becomes 

(~l) Sl +-+ s k aE sginfflT^j 1 ,^ 2 ) ... IT 1 ^, irf 1 ). (9) 

Lemma 14 When £ is small, and for every 6i,... ,6k £ {0,1}, 

£ k a E = (-i) k (-l) Sl+ - +s ’'dz{ 1 ...dz S k k . (10) 

Proof: 

Let UjVj be an edge of the path 7 ? where Uj precedes Vj. We can write 

Uj vj = e£(uj Vj) e ldj , ( 11 ) 

where £{ujVj) is the length of the edge UjVj in G, and 0j is the direction from Uj to 

Vj. Let us first consider the case of an edge UjVj oriented in the direction of the path, 
that is the dual edge w'jb'j of UjVj has its black vertex on the left of 7 1. By definition of 
the Dirac operator, we have K(w'j,b'j) = £(ujVj)e ld ie 1 2 . Next we consider the case of 
an edge UjVj oriented in the opposite direction, that is its dual edge Wjbj has its black 
vertex on the right of 7 Again, using the definition of the Dirac operator, we obtain 
K(wj,bj) = £(ujVj)e ld ie~ l 2 . We summarize the two cases by the following equation 

K(wf,bf) = (-i)(-l) s H(u jVj )e^. (12) 

When e is small we replace UjVj by dzj j . Thus combining equations (11) and (12) we 

obtain equation ( 10 ). □ 
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Lemma 15 When £ is small and up to a term of order 0(e), equation (9) equals 


( ~i) k 

(2n) k 


E 


[/ 




, s 2 b s l ( 0)1 


•••[/ ...dz s k , 


( 13 ) 


where F 0 (z,w ) = -, F\(z,w) = Fq(z,w), and the functions f w b are defined in 

z — w 

Section 5.2. 


Proof: 

Let us drop the superscripts Si- Plugging relation (10) in (9), we obtain 


(9) = (—i) k sgnaK 1 (bi, w 2 ) ■ ■ ■ K 1 (b k ,wi)^ u dz 1 ... dz k . 


| h - W j 
£—>0 £ 


Moreover, for every i ^ j, lim —-— = oo, so that by Theorem 11 we have 


K (hi, Wj) = —( F 0 (bi,Wj) + f Wj (0)Fi(bi, Wj )) + 0(e ). 


(14) 


(15) 


Equation (13) is then (14) where the elements K 1 (bi,Wj) have been replaced by (15) 
and expanded out. □ 


In what follows, all that we say is true whether the edge Wj J if I has its black vertex 
on the right or on the left of the path 7 ?, that is whether 5j = 0 or 1 . So to simplify 

notations, let us write {tj} instead of {5j G {0,1},^}, hence {tj} is the set of indices 
which run along the path 7 J. Keeping in mind that our aim is to take the limit as 
e —> 0, we replace the vertices bj and Wj in the argument of the function F ej by one 
common vertex denoted by Zj. Define 

H(e i,...,£fc)= 

= y] [fw2t 2 ^1*1 (®)] 1 • ' • [fwiti kfctfc k ^ £ 1 ’ Z 2t 2 ) • • • J z lt\ )dzit 1 • • • dzkt k . 


Lemma 16 


1. If (ei,... ,e k ) = (0,... ,0), then 


lim H( 0,.. 


Fq(zi, z 2 ) ■ ■ • F 0 (z k , Z!)dz 1 ... dz k . 


'71 


'Ik 


2. If (e i,...,£ fc ) = (1,...,1), then 

lim H(l,..., 1) = (-l) fc [ ... [ F 0 (zi,z 2 ) ■ ■ ■ F 0 (z k ,zi)dzi... 


dz k . 
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3 . Assume there exists i / j G { 1 ,..., k} such that £i = 0 , Ej = 1 , then 

lim \H(ei, .. .,£ k )\ = 0. 

£—>0 

Proof: 

Here are some preliminary notations. Dropping the second subscript, we consider an 
edge UjVj of one of the paths 7 ?, where Uj precedes Vj. Let us denote by Wjbj the dual 
edge of the edge UjVj , and let ee*“'' = v 3 — Wj, ee ■ = bj — v 3 . With these notations, 

we have dzj = £{e lolj — e* /3; ') (see Figure 4). 



Figure 4: Notations. 

Moreover, define 

J(£ 1, • • • ,£fe) = LWi (°)] £1 • • • [fwib k (Q)Y k dzi . ,.dz k . 

Proof of 1 . 

J(0,..., 0) = dzi... dz k , 

so that 

iL( 0 ,..., 0 )= ^ ...dz ktk . 

tl 

Since the paths 7 ^ are disjoint, the function Fq(zi, Z 2 ) ■ ■ ■ Fo(z k , zi) is integrable, and 
taking the limit as e —> 0 , we obtain 1 . 

Proof of 2 . 

.7(1,..., 1 ) ./ui2&i ( 0 ) • • • fwibk ( 0 ) dz\ ... dz k . 

Fix a vertex u of G*. Then, by definition of the function f wv , 

J( 1, •••,!) = fw 2 v(0)fvb 1 (0) ■■■ f WlV (0)f vbk (0)dzi...dz k , 

= fu>lV (0)f vbl (0).. ■ fw k v {0)fvb k (0) dzi ... dz k , 
fw\b\ ( 0 ) • • • fw k b k ( 0 ) dz\ ... dz k . 
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For every j, we have fwjbji 0) = e *(&+ a J'). Moreover, recall that dzj = e{e ia i — e* /3) ), 
so that — = e- l (Pi +a ^dzj, and we deduce 

= (-l) fc dzi ...dz fc . 

This implies, 

»(!.•• ■.!) = (-!)* £ Fo(z ltl , Z2t 2 ) ■ ■ ■ F 0 (z ktk , z ltl )dz Ul ... dz ktk . 

tl ?***j tk 

Taking the limit as e —> 0, we obtain 2. 


Proof of 3. 

Consider 0 < l < k, and assume E\ = ... = = 0, £(_ = ... = £ k = 1. Let us prove 

that lim e _>o |fL(0,... , 0,1,..., 1)| =0. Note that up to a permutation of indices, the 
argument is the same for the other cases. 

j(o,...,o,i,...,i) = f we+ Mo) ■ ■ ■ f wibk (0) dz\ ... dz k . 

As above, let v be a vertex of G*. Then, 


J(0,... ,0,1,... ,1) = (fvb e (0)dze)(f WlV (0)dzi)dz 2 .. .dzg-idzg+i . ..dz k . 
Introducing the following notation 


Hi = F(z 2 t 2 ,z 3t3 )... F{z t . 2 t t . 2 , ze-u^)F(z e+ltt+1 , z e+2tt+a ) • • • F(z k — ltk — 1 5 %ktk ) 5 
H 2 = F(zi tl ,Z 2 t 2 )F(z k t k ,z Ul ), 

H 3 = F(ze-i tt _ i: ze te )F(ze te ,ze + i tt+1 ), 

we obtain H(0,... ,0,1,..., 1) = 

= H 1 \ s ^H2fm 1 H v{0)dz ltl \ 53 H 3f vbete (0)dz Ue \ dz2t 2 .. .dz e _ ltt _ 1 dz e+ltt+1 ...dz ktk . 

t 2 ,...,t e ,...,t k \ ( i / \ *t ) 

Let us prove 

= °( e )- ( 16 ) 

te 

Dropping the second subscript, let ui,v\ = u 2 , v 2 = u 3 ,... ,v m -i = u rn . v rn be the 
edge-path y|. Denote by £ the quantity f m tl (0), then 

f vb . (0) = £e i(/3l " Ql) ... e~ iaj . 

Since dzj = £(e iaj — ), we obtain 


fvbj (0 )dzj 






JiPj-i-aj-!) 


_ _ gKdj-i-oij-i) e dPj-Qj)^ 
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hence 


fvbet e (fl)d z et e 

t g 


^ ^ fvbj (fydZj , 

J = 1 

m 

e £[l __ e *(/3l-“l) _|_ ^ e *(/3i-“i) _ _ _ e i(/3j-l-aj-l)) _ ( e *(/3i-“i) _ _ _ e i 0 3 j-i ~ a j-l) e KPj ~ a j 

i =2 

e£[l — (e*^ 1_ai ) ... e'Wm-am)], (telescopic sum). 


We deduce | f v b H (0)dzgt t \ < 2e, and (16) is proved. 

In a similar way we prove / Wl .u(0)dziti = 0(e). 

Using Taylor expansion in e for H 2 and H.$, we deduce that (X) tl H 2 f Wltl v(^)dz\t 1 ) and 
(j>2 tf (0)dz( t ^ are 0(e). Since the function H\ is integrable, we conclude that 

17(0,..., 0,1,..., 1) is 0(e 2 ) and so 3. is proved. □ 

Rewriting the second subscript, and summing equation (9) over the paths yf,..., , 

we obtain (by Lemmas 15 and 16): 

lim sgiid ^ E (-l) 8l+ - +Sk a E K-\b u{l , w 2tS ,)... K^fib^) = 
5i,...,4e{o,i} M Jk 

U 1 ‘“ u k 

= J2^ sgna ^ ( _1 ) £fc (/ "'/ F o(^i,4)•• ■ F o(4,4) d 4 ... dz £j , 

(17) 

where z® = Zi and z\ = Zi. When cr is a product of disjoint cycles, we can treat each 
cycle separately and the result is the product of terms like (17). Thus when we sum 
over all permutations with no fixed points, we obtain equation (5) of Proposition 12. □ 


Proposition 17 

- When k is odd, limE[(/i(uf) — h(u \))... (h(v f) — h(u%))\ = 0. 

6 —>0 

- When k is even, lim E[(/i(uf) — h(u\))... (h(v £ k ) — h(u £ k ))\ = 
k/2 

^ ' S(^t( 1) j W( 1) > u t(2) 1 v r(2) ) • • • d(^T(k— 1) 5 W(fc— 1)> ^r(fc) > W(fe)) j 
"rSTfc 

where g(u,v,u',v') = g(v, v') +g(u, u') — g(v, u') — g(u, v 1 ), g is the Green function 
of the plane, and T k is the set of all (k — 1)!! pairings of {1,..., k}. 



Proof: 

Let us cite the following lemma from [16]. 
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Lemma 18 [16] Let M = (my) be the k x k matrix defined by mu = 0, and mtj = 
t , , when i j. Then when k is odd, det M = 0, and when k is even 


det M = - -- 5 —— - 

rgTfe ( x t(1) “ x t( 2)) • • • ( x r(k- 1) 


x r(k )) 


Combining Proposition 12 and Lemma 18, when k = 2, we obtain 


limE[(/t(^) - h{u\)){h{v |) - h{u e 2 ))\ = 

£—>0 

= -T-o\ I [ 7-“— r^dzidz 2 + f [ jz — 1 _ dzidz 2 ) 

4 ^ yJ'Yi J 72 \ z i z 2 ) Jn J12 (^i ^2) / 

1 (fl - ^ 2 )(^1 - U 2 ) 

2ir 2 ° S (di - u 2 ){u\ - u 2 ) 

= -fl r (Ml,Wl,W2,'V2)- 
7 r 


The case of a general even k is an easy but notationally cumbersome extension of the 
case k = 2. □ 


5.4 Proof of Theorem 1 
Proposition 19 

limE [(H s <p) k ] = ±E[(F'p) k ] = 

£—>0 7r' 1. 

Proof: 

The second equality is just the /c th moment of a mean 0, variance ip), Gaussian 

variable. So let us prove equality between the first and the last term. 

Consider u\,..., distinct points of K 2 , and for every j, let vL- be a vertex of G £ lying 
within O(e) of Uj . Define 

H ij<P= e 2 a(v*)<p(v £ )(h(v £ ) - h(u £ j)) = ^ e 2 a(v*)ip(v £ )(h(v £ ) - h(u £ )), 
v e £ G e v e £K e 

where K £ = G £ n K , and K = supp(</?), then since we sum over a finite number of 
vertices, 


0 when k is odd , 

(k — 1)!! }, G(tp,tp) k / 2 when k is even. 


E[ H ^ip... H^ k p\ = E J2 - h{u\))... ]T e2a (. v k)‘P(. v k)( h ( v k) - h(u%)) 

vfeK £ v % eK£ 

= J2 { £2 ) ka ( v i)--- a { v t)‘P( v i)--- < p( v k) M K h ( v i)-h(u{))...(h(vl)- h {ul))]. (19) 

vf £K £ v?eK £ 
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Lemma 20 As e —> 0, the Riemann sum (19) converges to 


<p(vi)... <p(v k ) 


lim lE[(/i(vf) — h(u \))... (h(v |) — h{u k ))]dv\ ... dv k , 
£—>0 


where lim^oE[(/i(vf) — h(u \))... (h(v k ) — h(u k ))\ is given by Proposition 17. 


Proof: 

In what follows, all that we say is true whether Sj = 0 or 1, so to simplify notations, 
as before, let us write {tj} instead of {5j € { 0 , l},t^}, hence {tj} is the set of indices 
which run along the path 7 ?. Combining equations ( 6 ), ( 8 ) and (10) yields, 

E[h(vf)-h(ul)}...[h(vl)-h(u%)} = 

= (-*) fe s S naK l {bit 1 , w CT(i)t 0 . (1 )) • • • K ~ 1 (b ktk , wV(fc) t(r(fc) )^dzi tl ■ ■ ■ dzkt k ■ (20) 

ti.-.tko-gs; 


It suffices to consider the case where a is a fc-cycle, other cases are treated similarly. 
Indices are denoted cyclically (i.e. k + 1 = 1). There is a singularity in (20) as soon as 
Vj = Vj + i for some indices j. Hence, we need to prove that for e small enough, 

■ ■ ■J2( £2 ' )ka ( v i') ■ ■ ■ a ( v k)M v i) ■ ■ ■ <p( v i)\ l K_1 ( 6l *i > W2t 2 ). ■ ■ K _1 (6 fctfe , wi tl )\-^dzi tl . .. dz ktk , 

vf v% £ 

( 21 ) 


is o(l), when the sum is over vertices v \,..., v £ k that satisfy 

\vl ~v £ 2 \<6,...,\v £ m - v £ m+1 \ < 5, \v £ m+1 - v e m+2 \ >6,..., \v%- uf| > S, 

for some 1 < m < k — 1. Let 0 < /? < 1, up to a renaming of indices, this amounts to 
considering vertices vf,...,v k in ©1 n ©2 D © 3 , where 

©i = {v{,...,v £ m \ {1 < i < m, 0 < |< -< +1 | < £?}, 

@2 = {v £ m+1 , ■ • •, v £ n | m + 1 < j < n, s 0 < \v £ - v £ j+1 \ < 5}, 

0 3 = {v £ n+1 ,...,v £ k \n + l<£<k,\vf-v £ e+1 \> 6}, 

for some \ <m<n<k — 1. 

Since u \,..., u k are distinct vertices of M 2 , and since equation (21) does not depend on 
the path 7 ? from u £ j to v £ , let us choose the paths 7 f, ■ ■ ■ , 7 ^ as follows. Note that it 
suffices to consider the part of the path 7 J where the vertices bjt^ and '/177 are within 
distance 5 from v £ . Let us write \tj\ < <5 to denote indices tj which refer to vertices of 
7 j that are at distance at most 6 from v £ . Take 7 ? to approximate a straight line within 
O(e), from Vj to Vj + 5. Moreover, ask that if one continues the lines of and 7 £ +1 
away from uj and Uj +1 , they intersect and form an angle 6j. Let us use the definition 
of the paths 7 ?, and consider the three following cases. Whenever it is not confusing, 
we shall drop the second subscript. C denotes a generic constant, A ~ B means A and 
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B are of the same order. 


• v\,...,v £ m G 0i. 

By Remark 21 below, we have \K _1 (&*, it>j+i)| < C. This implies, 

1 / Q \ ^ 

/* ] \K (blti y w 2 t 2 ) ■ ■ ■ K (bmt rn yWm.+ lt rn+ i)\—^dzlt 1 ---dz m t rn < \~r~J 

* v n+V', v k e 0 3- 

By definition of the paths 7 |, \bg — wg+i\ > 6, so that lirrn^o &<! = oo. Using 

Theorem 11 yields ||/\' '(fy.«y..i)| = O (| 6<! _* £+1 |) < f • Hence, 

1 (c\ k ~ n 

^ ^ , ' w; n+2t n _|_ 2 ) • • • K (Pkt k ? ^lti )| fc _ n ^- 2 n+lt n _|_i • • • dZkt k < { — ) 

l*»+il<«.->l*fcl<* £ v y 

• G 0 2 . 

Let < L < 5 , and define the annulus, A(vj,L) = {u G G £ |L < \v — Vj\ < L + e}. 

By definition of the paths 7 ?, if v £ G A(vj +1 ,L), we have lim £ _,o —— 'f J+1 = 00 (since 
(3 < 1). Using Theorem 11 yields 

±\K~ 1 (b j ,w j+1 )\ =o( T , ---^ < —-^-rr-r = — C 


\bj-w j+1 \J min {w , +1 e7? } \bj-w j+l \ Xj + CL 


sin Oj , 


where Xj is the distance from v £ to bj. Let us replace sin 6j by C. Hence, if for 
m + 1 < j < n, Vj G A(vj + i, Lj), we obtain 

|L (fom+lt• • ■ LL {pnt n , m _j_i ■ • • dZnt n 7 


E 




E 


c 


c 




G CL m _|_l Xnt n G C L'i 


" • • • d,Xnt n 5 


5 G CL ; 


v CU 

j=m + l J 


c ](J log 


Let S be the sum (21) over vertices uf,...,u| G 0; n ©2 D © 3 . Denote by M = 
sup„eR 2 |¥’('y)|. Then, E < EiE 2 E 3 , where 


Hi = CAL" 


E £2 


E £ 2(fc - B) |^« + i)|... \<p{v%)\ 
••® c ©2 V 7 

5 S 

E -E 


"n + l--"fce 0 2 
i 

CM n ~ m ' 


_n — m r t 1 I ^ CLm-\-1 

£ J^m +1 • • • En log 


L m+ 1—0 L n —Q 


Lm +1 


... log I - + ] dL m+ i ... d,Ln 
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Moreover 


Si < 
“2 < 

S 3 ~ 


CM m 6 m £ m W- 1 \ 

E E 2(l "” ) k«+i)i'"iv(»i)i < v‘- n 

CM n - m (<51og5) n - m . 




Hence, E < CM k 5 2n k £ m ( 2 P B. Let us take (3 = 2/3 < 1, then rn(2fj — 1 ) > 0. If 

k-2n-\-l 

2n — k > 1, then E = o(l). If 2n — k < 0, take e < (W 2 * 3-1 ), and 5 = o(l). 


Remark 21 Let G* be a bipartite isoradial graph, and let K~ l be the corresponding 
inverse Dirac operator, then for every black vertex b and every white vertex w of G*, 
we have 

\K-\b,w)\ < C 

for some constant C which only depends on the graph G*. 


Proof: 

By theorem 4.2 of [18], K~ l is given by 

IF -1 (6, in) = J f W b(z) log zdz, 

where C is a closed contour surrounding cclw the part of the circle {e id \6 E [ 6 $ — it + 
A, #0 + 7 T — A]}, which contains all the poles of f w b, and with the origin in its exterior. 
Without loss of generality suppose do = 0. As in [18], let us homotope the curve C to 
the curve from —00 to the origin and back to —00 along the two sides of the negative 
real axis. On the two sides of this ray, logz differs by 2iri, hence 

K -\ b , w) = - / ^ Mt ) dt, where = ( , _ e .e, )(t _ ^ II 

Refer to [18] for the choice of path from w to b, that is for the definition of the angles 
0 \, 62 , otj.fjy These angles have the property that for all j, cos aj < cos (3j, so 


For every e ld G {e ld \ 6 G [0q — n + A, 0q + ir — A]}, and for every t < 0, we have 
1 1 — e l6 \ 2 > 1 1 — e l ^ +A ^\ 2 = 1 1 + e* A | 2 . Moreover for every t G M, we have 1 1 + e lA \ 2 > 
sin 2 A, and 1 1 + e* A | 2 > (t + l) 2 , thus 


[ \fwb(t)\dt < 

J — OO 



1 

¥TW 


dt T 



1 


sin 2 A 


dt 


1 + 


2 

sin 2 A 
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Hence |K 1 (b,w)\ < C, where C = ± (l + -4^-) 


We end the proof of Proposition 19 with the following 

Lemma 22 

1 . f ... [ (p(vi) .. .tp(v k ) limE[(/i(uf) - h{u\)) ... {h(v £ k ) - h(u £ k ))]dz 1 .. .dz k = 

0 when k is odd , 

(k - 1 )!!-t 7 ~G(<p, ip ) k/2 when k is even. 

2 . lim^oE [H^cp ...H^}= lim £ ^ 0 K[(H^) k ]. 

Proof: 

1. is deduced from the formula of Proposition 17, and from the fact that <p is a mean 0 
function. 2 . is a consequence of the fact that cp is a mean 0 function, and of estimates 
of the kind of those of Lemma 20. □ 

□ 


5.5 Remark 

Note that the double periodicity assumption of the graph G* is only required in Lemma 
13, where we implicitly use the expression of Theorem 3 for the Gibbs measure /i, as a 
function of the Dirac operator K, and its inverse K~ 1 . By Remark 4, this assumption 
can be released for the measure fi, in the case where the graph G is a lozenge-with- 
diagonals tiling. Hence, Theorem 1 remains valid when the graph G is any lozenge- 
with-diagonals tiling of the plane, periodic or not. 

6 Proof of Corollary 2 

We place ourselves in the context of Corollary 2: Q is the set of triangular quadri- 
tilings, and assume quadri-tiles are assigned the critical weight function. Recall the 
following notations: T £ is the equilateral triangular lattice whose edge-lengths have 
been multiplied by e, P is the Gibbs measure on Q of Section 2.4; for i = 1,2, Hfip = 
g2 Thv&v (T e ) ^r i P(. v )hi( v )i an d Tf , F‘2 are Gaussian free fields of the plane. 

In order to prove weak convergence in distribution of the height functions h\ and h ? 2 to 
two independent Gaussian free fields F\ and T 2 , it suffices to show that, V 99 £ C^(M 2 ): 

hm E[(HiV) fc (H 2 4n = E[(Ti^) fc ]E[(F 2 ^) m ]. 

£—>0 



□ 

□ 
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The key point is to obtain the analog of the moment formula of Proposition 12. The 
rest of the proof goes through in the same way, and since notations are quite heavy, we 
do not repeat it here. 


The idea to obtain the moment formula is the following. Recall that triangular quadri- 
tilings correspond to two superposed dimer models, the first on lozenge-with-diagonals 
tilings of C and the second on the equilateral triangular lattice T. Recall also that both 
lozenge-with-diagonals tilings and T are isoradial graphs. Hence, we start by applying 
Proposition 12 in the case where the graph G is a lozenge-with-diagonals tiling L 6 C. 
Then in Lemma 24, we prove some uniformity of convergence for every L £ £, and 
we conclude the proof by using Proposition 12 in the case where the graph G is the 
equilateral triangular lattice T. 

Let ui, ..., u k , v \,..., v k , ui,.. *, u m , Vi,..., v. m be distinct points of M 2 , and let 71 ,..., 7 ^, 
7 i,... , r / Tn be pairwise disjoint paths such that 7 j (resp. 7 ') runs from uj to Vj (resp. 
from u j to Vj). Define 


( 


\ 


G{u\,Vl, .. .,u k ,v k ) 


(~if 

( 27 r) fc 


E (-P' 




det 

i,fc] 



V ^ J 

Similarly, define (?(ui, vi,..., u m , v m ). Let Uj,Vj, Uj,v^ be vertices of T £ lying within 
0(e) of Uj,Vj,Uj,\/j respectively. 


Lemma 23 


limE[(/ii(uf) - MO).. . (/ii(vfc) - MOXMO - h 2 ( uf)).. . (MO) - MO))] = 
£—>0 

= G(ui,Vl, . . . ,Uk,Vk)G( Ul, Vl , . . . , U777,, V777,). 


Proof: 

By definition of the measure P, we have: 

H(hi{v f) - MO)... (MO - MOXMO - Muf)) • • • (MO) - Mi4))] = 

= E M - M«D) ■ ■ ■ (M<4) - M«!))(M v !) - Mu?)) • • • (M v m) - M^m))] V(l*), 

L*€.M(T*) 

E (M*!) - W)) • •' ) - M« 5 )) • • • (hi(v%) - hiK))]V(L*). 

L*gM(T*) 

Using Section 5.5, we can use Proposition 12 in the case where G is any lozenge-with- 
diagonals tiling L E C (periodic or not), hence for every L E £, we have: 


iimE z,[(/ti(uf) - MO) • • • (MO - MO)] = G(ui,v 1 ,.. .,u k ,v k ). (22) 

£—>0 ^ 

Note that the right hand side is independent of L, hence to obtain Lemma 23, we 
need to prove that convergence in (22) is uniform in L, see Lemma 24 below. Indeed, 
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assuming this is the case, for z small, we can write: 

Y (M v l) - Muf )) • • • (M v m) - h 2 (u S m ))m tl L [{hl(vl) - hl(uf)) . . . (fcl(«k) - ^l04))] V(L‘) = 
L*eA1(T») 

(M v l) - M u l)) • • • (M v m) - ft 2 (u^))(^(ui, Vl, Ufc) + O 0 M‘ T ( L*), 

L*6«(T*) 

= (P(ui,Wi,...,Ufc,Vfc)+ 0 (e:)) (M v !) - /i 2 ( u D) • • • (M v m) - h 2 (ii £ m ))dy T (L*), 

L*eA1(T*) 

= (S(«i, ui, • • •, u k , v k ) + 0 (e)) E^t[(/ i 2 (vf) - ft 2 (uf))... h 2 {y e m ) - ft 2 (u^)], 

= (G{ui,vi,.. .,u k ,v k ) + O(e)) (G(u 1 ,v 1 ,..., u m , v m ) + 0 (e)), 

where the last line is obtained by using Proposition 12 for the graph T. 

Lemma 24 When e is small, and for every lozenge-with-diagonals tiling L e £, 

- hiK))... OiOI) - M*4))] = • • -,Uk,Vk) + O(e), 

where 0(e) is independent of L. 

Proof: 

Let us look at the proof of Proposition 12 in the case of a lozenge-with-diagonals tiling 
L £ C, and denote by K k the Dirac operator indexed by vertices of L*. Then Lemma 
24 is proved if we show that O(z) in Lemma 15 is independent of L £ C. Looking at 
the proof of Lemma 15, we see that 0{e) comes from the error term in the asymptotic 
formula for the inverse Dirac operator of Theorem 11, [18]: 

Kl l (bi,Wj) = e wj) + f Wjbi (0) F i( b i , wj )) + 0 ( _ £ )) . 

In [18], the error term is computed explicitly. Looking at the explicit formula, and 
using the regularity of the graphs L, we show that (){ \ b f w \i ) = where C\ is 

independent of L. Moreover, by assumption the paths 71 ,... are disjoint, so that 
we define: 

C -2 = inf inf \b — w\ > 0 , 

* 7 £j {b£'/i,w£'yj} 

which is independent of L. Hence ()( \ b J w \ 2 ) = ■ □ 

□ 
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